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Abstract 

The decay process B — > D + D~ir° is an interesting channel for the investigation of CP 
violating effects in the b— sector. We write down a decay amplitude constrained by a low-energy 
theorem, which also includes the contribution of resonant S— and P— wave beauty and charmed 
mesons, and we determine the relevant matrix elements in the infinite heavy quark mass limit, 
assuming the factorization ansatz. We estimate the rate of the decay: B(B — ^ D + D~tt°) ~1x 
10~ 3 . We also analyze the time-independent and time-dependent differential decay distributions, 
concluding that a signal for this process should be observed at the B-factories. Finally, we give 
an estimate of the decay rate of the Cabibbo-favoured process B — ► D + D~K$- 
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1 Introduction 



Multibody hadronic B decays represent a large fraction of the inclusive nonleptonic rate [U, and 
therefore it is worth analyzing their phenomenological aspects, since they constitute accessible 
channels for the experimental investigations ]2|. In particular, some three-body neutral B hadronic 
decays deserve further attention, from both the theoretical and the experimental sides, since they 
have been recognized as an important source of information concerning CP violation in the beauty 
sector. This is the case of the decays S° — > ir + ir~ir° and B° —* ir + ir~ir°, which provide information, 
together with the neutral B decays into a pion pair, on the CKM angle a |, 1. This is also the 
case of the three-body decays i?° — > D + D~tt° and B° — > D + D~tt°, which have been identified as 
interesting channels to investigate the angle Q, in particular as far as the discrete ambiguity 
(3 —> 7j — (3 of the CP asymmetry in B — > J /ipKg is concerned. The removal of such ambiguities 
and, in general, the identification of possible constraints on the CKM angles are of prime interest, 
mainly in view of testing the Standard Model and probing the effects of new physics scenarios. 

S3 

The theoretical calculation of multibody hadronic B decays presents uncertainties, mainly re- 
lated to the long-distance QCD effects involved in these processes. Simplifying assumptions are 
usually adopted, such as, for example, the hypothesis of dominance of intermediate hadronic res- 
onances in the relevant amplitudes. In the case of pions in the final state, however, low-energy 
theorems can be employed to reduce the decay amplitude in the soft-pion limit q n — > 0; this allows, 
for example, to relate a three-body decay amplitude to a corresponding two-body one. If a narrow 
phase space is available around the point q n — > 0, an extrapolation can be done to estimate the 
multibody process. This program cannot be pursued for the B — > 3n decays, where high momen- 
tum pions are allowed in the final state. The situation presents less difficulties in the case of the 
decay B —* D + D~ir°, where a quite narrow phase space is available for the pion; therefore, the 
amplitude having the right behaviour for q w —* can be extrapolated, including the contribution 
of intermediate resonant states, to the full phase space. 

This is the aim of the present work. We shall write down an amplitude for B — > D + D~ir° and 
for the SU (3)-related process B — > D + D~ K$ having the soft pion limit required by current algebra 
and PCAC, and including a set of intermediate hadronic states. In this way, the amplitude can be 
reduced to a set of two-body hadronic matrix elements, which we shall evaluate by the factorization 
ansatz; the description will be simplified by observing that, in the infinite heavy quark (6, c) mass 
limit, the hadronic matrix elements involved in the calculation are related to a few universal (mass- 
independent) parameters. An interesting observation will be that the full amplitude can be derived 
from an effective Lagrangian, obeying chiral symmetry in the light meson sector and heavy quark 
symmetry in the heavy quark sector. The unknown parameters are the Isgur-Wise semileptonic 
form factors, the heavy meson leptonic constants and the effective couplings describing the QCD 
interactions of the heavy mesons with pions. 

The plan of the paper is as follows. In Section |2| we briefly review the kinematics of the 
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B — * D + D~7T° decay, and the relevance of this channel in the perspective of CP measurements. 
In Section ^ we derive the low-energy theorem for the nonleptonic B — ► D + D~ir° amplitude, 
together with the contribution of intermediate resonant states, and provide an evaluation of such 
an amplitude by using the factorization ansatz. In Section [| we discuss a derivation based on an 
effective heavy meson chiral Lagrangian. The numerical analysis is reported in Section ||, and a 
short discussion of the Cabibbo-favoured B — * D + D~ K$ decay concludes the presentation. 



2 Kinematics and (3 dependence 

We consider the process: 

B°(p)^D+(p + )D-(p„)Tr°(q) (1) 

and the analogous one for the B° meson. Neglecting penguin contributions, these decays are 
governed by the weak Hamiltonian 

HW = 7§ y ^(ci + ^)67 m (1-75)cc 7m (1-75M +h.c. ; (2) 

Gf is the Fermi constant, V{j are CKM matrix elements and ci^ are short-distance coefficients, 
with N c the number of colours. 

The neglect of gluon and electroweak penguin operators, that in principle contribute to process 
(|l]), cannot be justified a priori. However, considering that the corresponding short-distance coeffi- 
cients are rather small, and that process ([!]) is colour-allowed, we may expect the dominance of the 
"tree-diagram" operator in Eq.(^) to be a good approximation in this case. Qualitative estimates 
based on the factorization ansatz suggest for the two-body b — ► ccd decays that corrections from 
the penguin contributions could be of the order of few percents [J5J. We recall that, due to the 
common AI = 1/2 character of both tree and penguin operators for the transition b — ► ccd, the 
corresponding amplitudes cannot be separately determined by the isospin analysis §§. 



Following the notations of Ref . 0] , we define the Dalitz plot variables of the decay (|T|) : 

s+ = {p+ + q) 2 

s_ = {p-+qf (3) 
S = (p+ + V- ) 2 = m B + 2m Z) + m l - s + - s - 

In terms of the heavy meson four-velocities 

pt 1 = niBV^ , pi = mr)V+ , p^_ = mj)t)f , (4) 

we also introduce a set of invariant variables, suitable for the application of the heavy quark effective 
theory (HQET) §§ to our problem: 

s+-m 2 D - ml 
q + 2m D 
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q ■ V- 



i i 
m D — m%. 



2m d 



s+ + S- — 2m 2 



v + ■ V- 
v ■ v + 

V ■ V- 



2niB 
m 2 B + m 2 — s + 

2m% 
m 2 B + m 2 D — S- 

2m£,mB 
m 2 B + m 2 D - s + 



(5) 



In the plane (s_, s+) the allowed kinematical region is bounded by the curves 



sT n = (E 2 + E 3 y-[JE 2 -m 2 n + JE 2 -m 2 D (6) 



(E 2 + E 3 ) 2 - JE 2 -ml- Je 2 - m 2 D (7) 



^ s~—m 2 n + m 2 , ^ — s_ — m 2 n + m 2 , , . , 9 . , 9 

where E 2 = - - and E 3 = — , with (mo + < S- < (tub — mo) ■ 

-\ ~ -\ ~ 

The kinematical region is symmetric under the exchange s_ «-> s+; CP eigenstates correspond to 

the line s + = s_. 

The role of three-body decays in accessing the weak angle (3 has been discussed in Ref. §, and 
we repeat here the basic points and the notations of that analysis. Since no weak phases appear 
in the effective Hamiltonian (|2|) in the Wolfenstein parametrization, the only relevant weak phase 
in the B° (and B ) decays of the kind (|l|) is the phase (3 of the mixing B° — B . Denoting, as m 
Ref. by A(s+, s_) and A(s + , s_) the amplitude for the decay into D + D~ir° of the B° and B , 
respectively, the time-dependent decay probability of a state identified as a B° at t = is given by: 

p-rt 

-» D + D-it°)\ 2 = — -[G (s + ,s_) + G c (s+,s_)cos(Amt) 

- G s (s+,s_)sin(Amt) , (8) 

where 

G (a+,O = Ws +!S -)| 2 + p( S+!S _)| 2 
G c ( s+ , s _) = \A(s + ,s.)\ 2 -\A(s + ,s.)\ 2 

G s (s + , s_) = -2 sin(2/3)i?e(.A* (s + , s_)^(s + , + 2 cos(2/3)Im(.4* (s+, s_)^4(s+, «_)) 

= -2sm(2p)ReG s (s + ,S-) + 2cos(2f3)ImG s (s + ,s^) . (9) 

For the analogous decay of the B° one has 

_n e~ Tt r 

(t) -»■ Z? + Zr^°)| 2 = — [G (s-,s+) -G c (s_,s+)cos(Amt) 

+ G s (s_,s+)sin(Amt)j . (10) 
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Assuming that no direct CP violation occurs, consistently with the neglect of penguin operators, 
the condition A(s + , s_) = *4(s_, s + ) is verified; in this case, the time-independent term Go(s+, s_) 
in Eqs.(p|)-(p!C|) is symmetric in s+ <-> s_, while the coefficient G c (s+, s~) of the cos(Ami) term 
is antisymmetric. In principle, such contributions to the decay rate can be directely tested by 
symmetric or, respectively, antisymmetric integration over of the time-dependent Dalitz 

plot distribution of events. 

As far as the interference term sin(Amt) in (p|)-(|lO|) is concerned, G s (s+, s_) would be symmetric 
under s+ «-> s_ in the case of real A(s+, s_), and only the "indirect" CP violating part proportional 
to sin(2/3) would survive. In principle, also the effect of this term can be disentangled from the 
other ones by symmetric integration in s + and s_ of the experimental time-dependent Dalitz plot. 
In the general case, however, the amplitude j4(s+, s_) will have a non-trivial CP-conserving phase 
5(s+,s-), of strong interaction origin, such that both the CP-violating sin(2/3) term and the CP- 
conserving cos(2/3) term will contribute to Cr s (s+, s_). In particular, the role of the latter term 
was emphasized in Ref . Q , following the discussion of the decay process B —> 3tt of Ref . ||] , as a 
possible resolution of the discrete ambiguity (3 — * ^ — (3 implicit in the experimental determination 
of sin(2/3) from, e.g., the time-dependent CP asymmetry of the process B — * J/tfjKs- One can 
easily see that, in the hypothesis of no direct CP-violation, the sin(2/3) component of G s (s+,s_) 
should be symmetric under s+ <-> s_ as being proportional to cos(<5(s+, s_) — 5(s_,s+)), whereas 
the cos(2/3) term should be antisymmetric as being proportional to sin(<5(s+, s_) — 5(s_ , s+)). Thus, 
under such assumption they could be disentangled by symmetric and, respectively, antisymmetric 
integration of the sin(Amt) component of the Dalitz plot distribution. 

According to the above considerations, Eqs.(|8]) and ( |TCj| ) imply that a time-dependent analysis 

of the neutral B decay to D + D~tt° gives access to cos(2/3) if the product A* A has a non-vanishing 

imaginary part. Clearly, the required CP-conserving strong phase between A(s+, s_) and A(s + , s_) 

must have a non-trivial dependence on s+ and s_. Following Ref. [Q], we assume the variation 

of such strong interaction phase over the Dalitz plot to be entirely determined by a set of excited 

B* and D* resonance contributions, parameterized by Breit-Wigner poles in the relevant channels. 

In addition, however, considering the rather low energies (on the heavy quark mass scale) allowed 

to the pion in the considered decay, we constrain such polar expression for the decay amplitudes 

to obey the low energy theorem resulting from chiral symmetry. In the appropriate limit, q — > 0, 

the amplitude reduces to the continuum "contact" term determined by the general (and model 

independent) current algebra procedure. Unavoidably, the assumed resonance behaviour in s + and 

s_ , as well as the factorization approximation for the relevant two-body matrix elements, introduce 

some amount of model dependence that is difficult to reliably assess on purely theoretical grounds. 

On the other hand, the experimental study of the Dalitz plot distribution of events should allow 

to test the phenomenological validity of the model, and in particular to evidence non-resonant 

contributions to the strong phase variation if they turned out to be large.Q In any case, we shall 

"The sources of systematic uncertainties implicit in the assumed resonance parameterization of the strong phase 
behavior have been discussed in detail for the three-body B — » pn — > 3-7T decay in pp. 
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estimate the theoretical uncertainties of our approach by considering two different extrapolations 
for the Breit-Wigner poles. They will be discussed in the next Sections. 



3 Low energy theorem and polar contributions 



In order to derive a low energy theorem for the amplitude (|l|) we consider the Ward identity [ 10 1 



A 



< D + {p + )D-{p^(q)\H w \B\p) > 

-ill < D + (p + )D-(p^\[Fi,H w ]\B°(p) > +lim q ^ 

Jtt 

Ab + Ar. 



il r q X M x - A B 

Jtt 



In Eq. flllD, M\ is 

M x = i J <fxe^ x < D + (p + )D-(p„)\9(x )[A x (x),H w (0)]\B°(p) > , 



fill 



(12) 



with A\ = (u^fxj^u—d'jx'y^d) /2 and F§ the corresponding axial charge; H\y is the weak Hamiltonian 
in Eq. (||) and f n = 132 MeV. The polar contributions As and Ar in Eq. (|Tl|) are depicted in 
fig.l. The first set of contributions Ar includes those intermediate states which become degenerate 




D D 



(a) (b) 

Figure 1: Polar diagrams contributing to Ar and Ar. The dot corresponds to a strong vertex, 

the box to a weak vertex 

in mass with the initial B or the final D states in the HQET: the J p = 1~ B* (fig. la) and D* 
(fig. lb) mesons. The second set Ar denotes contributions from excited beauty and charm mesons, 
corresponding to P- waves in the constituent quark model: Bq, B% and Do, with J p = (0 + , 2 + ), 
respectively. Clearly, this is a simplification, since in principle the contribution of other intermediate 
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resonances can be considered, such as, e.g., B — > ip^n — > D + D~ir°. In the case of ip' , which 
should be the most important one, an estimate of this colour-suppressed process in the factorization 
approximation gives a negligible contribution with respect to the other ones considered here. 



Since H\y has a (V — A) x (V — A) structure, for the equal-time commutator in Eq. ( [LTD the 
equality [F^,H\y] = — [F 3 ,H\y] holds, with F 3 the isotopic spin operator. Then, using F 3 \B° >= 



\B° > and F 3 \D+D- > 



(J=l,0;/ 3 =0) = 



0, the equal-time commutator becomes: 



Ai 



-i^L < D+(p + )D-(p.)\H w \B°(p) > 



(13) 



The separation indicated in the second and the third terms on the right hand side of Eq. (|11| ) is 
done in order to avoid the ambiguity which arises in taking the mass degeneracy limit first and 



then the limit q — > or viceversa in iq x M\ or Ab, so that — i^p-q x M\ — Ab 

L J n 

limit when q — > 0. Ab and M\ (Born) (which alone is relevant for the above 
calculated, and clS ct result Eq. (|TT[) becomes: 



has a well defined 
imit) can be easily 



■B* 



r D* 



A(s + , s_) = Ai(s + , s_) + A 2 (s + ,s^)+A 4 (s 



+ A 5 {s + ,s.)+A 



R 



where 



A 2 (s+,s-) 

^4 j5 (s + ,S_) 



r 1 



gB*0BOn 



1 + 



m. 



m 



D 



1 + ■ 



so 



(P - ff)t 



m 



B* 



so 



-(p + ?)/ 



(14) 



(15) 



m 



D* 



D* 



{p± + q)ix + 



m 



<P±-q)v G» ± (16) 



D* 



with 



F»e, 



< D + D~\H W \B 

< D ± D*\H W \B° 



*o 



> 



> 



(17) 



e and rj being the B* and D* polarization vectors, respectively. It may be noted that the constant 
terms in Eqs. ( |l5| ) and (16) correspond to the limit q — » indicated in Eq. (O). Here, consistently 

2 2 

with the use of the infinite heavy quark mass limit, we have neglected terms of order 5b = B 2m B B 
and 5b in comparison with the heavy meson masses. 



With H\y given in Eq. (g) and using the factorization ansatz, the matrix elements ( |13|) and 

(|i~7|) can be evaluated in the heavy quark effective theory (HQET) in terms of the Isgur-Wise 

semileptonic form factor £(v ■ v'), where v and v' are the relevant four-velocities. As for the polar 

terms in Ar in Eq. (|j~4|), we only consider P-wave intermediate charm and beauty resonances. In 

this case, the matrix elements corresponding to the relevant amplitudes can be written in terms of 

G c 

the Isgur-Wise universal form factors 7i/ 2 (i' ■ v') and t 3 / 2 (v ■ v'). Defining K, = —=(c\ + -jrr)VcbV^, 
and parameterizing the effective strong couplings 

1 _ 2m B 

1 



9b*+b ^ 



9b*°b°-k° 



9d*-d°tt- 



V2 



9d*-d- 



9 



-9 



(18) 
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and the current-particle vacuum matrix elements 



< £r(p_)|cy(l - 7 5 )d|0 >= i- 



<D (k)\c^(l- 75 )d\0>-- 



—=pr 

F+ 

-i , 



(19) 



we can list the expressions of the various contributions introduced above. The equal time contri- 
bution is simply given by: 



-r < \ fCFJmB ^,m 2 B + m 2 D - s- .rm 2 B + ml - s. 
Ai{s + ,S-) = j=- — ) 



V2f n 



2m B m B 



2ni£) 



+ 



m 



D 



m 



B~ s + 



2m B 



(20) 



As for the polar terms, the contribution of the B* intermediate particle in fig. la is: 



\[2KFg Jmsm b* 



so 



A2(S+,S-) = " """ " V "2 g(- 

f w s - m%, 2m D m B * 



so 



2m£>mB* 



+ 



s s + - m 2 D - ml 



so 



(21) 



2m,B* 2ni£) 

(we neglect the B* width); on the other hand, the contribution of the Bq, J p = + state, whose 
width is Tfe, reads: 



A3(s+,s_) = 



K,FGb Btt(so) 



1 



:Tl/2(; 



so 



-)[(m Bo - m D )s - 2m 2 D m Bo ] ■ (22) 



«o ~ m B + irn B () ^b 2rn Dy /m Bo ' ± ^ y 2m Bo m D ' 
The contributions of the poles D*~ and D* + in fig. lb (Yd* is the D* width) are, respectively: 

s + + s_ — 2m 2 D 



A 4 {s + ,s-) = 



A 5 (s + ,s-) 



+ 



\[2KFg y/m~B~m£)m£)* ^m B + m 2 ^ — ^m B + rri£) 



m 2 D , + imo^D* 



2m B m£> 



m B m D 



s_ — m 2 D + m 2 m 2 B — m 2 D + s. 



2?7T.£j* 



(23) 



V2ICFg 



^ m 2 B -m 2 D + 



fw 4^/niB s+ — m 2 D * + irriD*TD* v 2m B 'mD* 
— (s+ + s- - 2m 2 D )(m 2 B - m 2 D - s + ) + (2m B m D » + m% - m 2 D + s+)(s_ - m 2 D - m 2 ) 

(24) 



m B 1 2 2 \/ 2 2 1 \ 

\m B -m D - s + ){m w - m D + s+) 



m D * 



The contributions of D and Dq in fig. lb can be written as 



Ae(s + ,s-) 
and 

A7(s+,s-) 



KG DoD „{s-)F + 



-di- 



ll 



2^/m Do m B m D s_ - m 2 Do + im Do T Do ^ y 2m B m D 
1CG DoDn (s + )F 1 



-){mB-m D )[(m B +m D ) 2 -s. 

(25) 



2^/m Du m B m D s+ - m 2 Do + im Do T Do 
m B {m 2 B - m 2 D - s+) 



Tl/2( 



2m.Bm Do 



)[m Do (m 2 B + m 2 D -s + ) 
(26) 
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respectively. In the previous equations, the following definition holds: 



Gd d-k(s) 



m Do m D s 



m 



D 



h 



2 m Do f n 
and an analogous expression is used for Gb Btt- 

Finally, we consider the contribution of the pole, whose width is Td 2 : 

/rnB~m D2 (m B +m D2 ) 



(27) 



—t i \ /R h ,s + + m%-m 

As(s+,s-) = OV6— r 3/2 ( 

U 2m B m D2 

{ 

1 

3 



s + + S- 



2m 



D 



+ m 



m 2 D2 + im D2 T D2 



B 



m 2 D s + 



2m b 



2mBm£) 2 



+ m 



B 



111 



D\2 



2mBm£) 2 

and the contribution of the B 2 ° intermediate state: 



m„ 



2m D2 
m 2 D + m 



2m D2 



1 1 1 
^) 2 ]} 



(28) 



Ag(s + ,S-) 



-ICFV6^t 3/2 (—^ 

f n 2m B m B2 



fm^m B2 {tub + m D2 ) 



so ~ m B2 + imB 2 TB 2 



m 



D 



mt 



so 



s + + s_ — 2m 



2ml 



2m£) 



2mBfnr> 



2m B2 



1 



Am B m B2 



+ s. 



2m 



D 



2m 



1 1 1 



j ( 29 ) 
2m B2 J J 

The D 2 ~ contribution vanishes in the factorization approximation. Notice that, for simplicity, we 
have assumed momentum-independent widths in the Breit-Wigner denominators. 

In the above equations, the usual definitions of the universal Isgur-Wise form factors have been 
used (see, e.g., the reviews [^, |ll|); as for the the effective coupling h! in the D 2 Dir and B 2 Bir 
vertices, it has been first investigated in [12] in the framework of HQET and we shall turn to this 
coupling in the next Section. 

Eqs.(^l|)-(p^) are obtained by considering the expressions for the effective strong vertices and 
the weak ones in the factorization approximation, and combining them to evaluate the diagrams 
in fig. la, b. This procedure presents some uncertainties, for example related to the relative signs 
between the various contributions. A method which allows to partially overcome such difficulties 
is based on the use of an effective heavy meson chiral Lagrangian, and the next Section is devoted 
to this approach. 



4 Evaluation by an effective chiral Lagrangian 

In order to determine an expression for the amplitude (Q) let us consider the effective Lagrangian 
PI 



C = igTr(HH^ 5 A^ + ihTr(HS^ 5 A„) + -r- ^{HT^^D^Ax + h 2 D x A fl }) + h.c. 
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that describes the strong interactions of pions and kaons with heavy mesons containing one heavy 
quark. This construction of the effective vertices follows the prescription of HQET, with the further 
constraints imposed by chiral symmetry. H, S and T represent heavy meson doublets corresponding 
to different values of the spin-parity sf of the light degrees of freedom of a qQ meson; the doublet H 
comprises the negative parity low lying states, viz. D, D* in the case of charm and B, B* for beauty; 
the multiplet S is characterized by sf = ^ + and comprises the positive parity (J p = + , 1 + ) low 
lying states, viz. Dq, D*' for charm and Bq, B*' for beauty; the multiplet T has sf = | + and 
comprises the positive parity (J p = 1 + , 2 + ) states: D\, D\ for charm and B\, BZ, for beauty. 



The fields H and S and T are 4x4 matrices containing annihilation operators. In the charm 

l- 

2 



sector, for the negative parity states sf = A these fields are given by 



H = [jy^f - D 1B ] , (31) 

and the conjugate field is given by H = jqH^jq. For positive parity sf = (\ + , § + ) states, the fields 
are defined by 

S = ^[Dfcf-ys ~ Do] (32) 



rpfl 



(33) 



In Eqs.(|31~|)-(|33|) v generically represents the heavy meson four- velocity, D*^, D, D] 1 ^ and Dq are 
annihilation operators normalized as follows: 

(0|Z>|cg(0-)) = ^fm ^ (34) 



(0|D*"|cg(l-)) = e^VMn , (35) 



and similar equations hold for the positive parity states (in Eqs. (p4|) and fl35|) Mh = = Mp* is 
the common mass in the H multiplet); the transversality conditions are v^D*^ = v^Df'^ = v^D*^ = 

The couplings HHtt, HStt and HTir of the heavy mesons with light pseudoscalar mesons are 
constructed through the axial vector current 

^ = l(^_^t) ; (36) 

where £ = exp(ill// 7r ), with II the familiar 3x3 SU(3) matrix describing the octet of light 
pseudoscalar mesons. As it is clear from Eq.(^), the interaction vertices HHir, HStt and HTtt are 
described in terms of the effective couplings g, h and hip. We shall quote in the next Section the 
numerical values for such parameters; here, we only notice that in our calculation the combination 

h' = ^ (37) 

(A x is a mass parameter) is needed, which can be determined from the experimental measurement 
of the Z?2 pionic transitions. 
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In terms of the heavy and light meson operators, an effective weak nonleptonic Lagrangian can 
be written as follows: 

C eff = ^ (ci + V cb V: d Tr[(L M + L'») (S) (J» + ^) (c6) ] • (38) 
The effective currents , (the subscript c means anti-charm) are 

Lfa = l -F ^(1-75)%)^) (39) 



Lg, = i F+Tr^l - 75)%)£ f ) (40) 



with i* 1 and -F + already introduced in (|i~9|). The effective and currents describing the weak 
b —* c transition can be written in terms of universal form factors: 

= -^(^•^)Tr(%)K)7 /i (l-75)%)W) (41) 
J'fo = -Ty 2 (v-v')Tr(S {c) (v')r(l-^)H {b) (v)) (42) 

J "cb) = "V^-3/2 (^^)^(^T ( A C) (^)7 M (1-75)^)M) (43) 
where t>^, i/^ are the heavy meson four- velocities in the initial and final state. 



The Lagrangian (38) meets the following requirements: i) it allows b — > ccd transitions with B 



(or B*) in the initial state, and two charmed mesons (with sf = ^ , \ or | ) plus any number of 
pseudoscalar light mesons in the final state; ii) the resulting amplitude corresponds to the evaluation 
of the weak 4-quark effective nonleptonic Lagrangian in the factorization approximation; iii) it 
contains the minimum number of light meson field derivatives, consistently with general properties 
of chiral symmetry and the soft-pion limit procedure employed in the previous Section. 

The effective weak nonleptonic Lagrangian (^) , together with the strong interaction Lagrangian 
(p0|), allows to write down an expression for the set of amplitudes contributing to the transition 
(|l]). The equal-time contribution derived from Eqs,( |30"| , |3"8| ): 

Ai(s+, O = ^ F £{v ■ v + )(v + . v _+vv^) (44) 

exactly reproduces Eq . (pc|) taking into account the invariants in Eq. @. Together with this term, 
the set of polar contributions corresponding to Eqs.(pl|)-(p9|) can be written; the differences with 
respect to the expressions reported in the previous Section represent a set of finite mass corrections, 
that partially account for the theoretical uncertainties in the calculation. They correspond to 
different treatments of the Breit-Wigner forms. 

The B* and Bq (of width T^) pole contributions are given by 

—a , x JC^W^m D ~ £(v ■ v+) 

A 2 {s + ,s-) = Fg 



a/2/tt q-v + {m B * - tub) 

( v ■ t>_ (q ■ v + — v ■ v + q • v) — (1 + v ■ v+)(q ■ V- — v ■ v _ q ■ v)) (45) 
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and 

K,^fm~^m D a , , q-v+ (-v+ • v- + v • v-) , . 

A 3 {s + ,s_) = j=- h F t 1/2 {vv + ) — , (46) 

V^/tt —q-v — Om + i-f 

respectively, with 5m the mass difference 5m = Mb — Mb = Md — Mjj. As for the D*~ and D* + 
contributions (Td* is the D* width), they read respectively: 

—t i v K,^m~^m D ~ q-v + +q-v-q-v-(v + -v-+vv-) 

A A (s + ,s_) = gF£{v-v + ) T - , (47) 

V^/tt q ■ v- - (mo* - mo) + 

—jr , x K,Jm^m D a £{v • f+) 



VZU q-v+ — (mo* - mo) + 

(1 + v ■ v + )(q ■ V- — v + ■ V- q ■ v+) — v + ■ w_ (q ■ v — v ■ v + q ■ v + )] . (48) 



The Dq and J p = + pole terms are 

a/2/tt q • v- — 5m + i 2 



a-) = ^ " h F+ Z(?> ■ v+) ~ r .r, ( 49 ) 



and 

■47(a+,a-)= at, /i J iT ri/ 2 («-v+) " rf- • ( 50 ) 

v^/tt g ■ d + - Am + if 

The D 2 + pole contribution reads: 

-r , x KJm~B~m D ^h'F 

As{s+,s-) = v t=- — p- t 3 / 2 (v«+) 

V2/7T c7 • - 0m 2 + z-^ 

/ V ■ V + + l r n , 91/ -v 

X ( 3 Kr- (?• U+) ](-«+ + V-V-) 

+ [q-v-(v-v + )(q-v + )][(v-v + + i)q-v--v + -v-(q-v + + q-v)]) , (51) 

with 5m2 = Mjj* — Mp = Mb* — Mb- Finally, the B 2 ° pole contribution is 

-r , x KJm~B~m D y/3h'F 

Aq[s+,s-) = 7 =- — t -t 3/2 {v-v + ) 

V 2/tt —q-v — 5m 2 + z-j 2 

/ U • f + + 1 r o / \2i/ \ 

x ( g [m ff -(g-i;) ](«+•«_ -«•«_) 

- [g • n + — (u • • u)]^ • u+ + 1)<? • v- — v ■ • v + + g • v)fj. (52) 

The expressions for the various terms Ai — Ag allow us to reconstruct the amplitude describing 
the process (|l]). The differences with respect to the results obtained from the amplitude derived 
in the previous Section represent a theoretical uncertainty associated to the writing of the polar 
contributions. 



5 Numerical analysis 

In order to estimate the rate of the decay (||) and the terms in Eq.(j9|), using the formulae in the 
previous Sections, we must rely on numerical values for the various hadronic parameters such as the 
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leptonic constants, the strong couplings and the semileptonic form factors. In some cases, experi- 
mental information can be used; theoretical methods can be adopted to determine the remaining 
quantities, and for these we shall mainly use the results of the QCD sum rule method. 



The strong coupling constants appearing in (^C]) have been evaluated by several authors |13|, 14 1; 
we assume here for g and h the values given in Ref. |]|: g = 0.35, h = -0.52. The combination 
h' = hl + h2 can be obtained from experimental data; from the full width of the .D 2 (2460) state, 
r 2 = 23 ± 5 MeV g], we get ti = 0.60 GeV" 1 , assuming that D\ Dtt saturates the hadronic D\ 
width. 

The constants F, F + do not depend on the heavy quark mass (modulo logarithmic corrections) 
and have been estimated by QCD sum rules ||, m, H|. We take the values: F = 0.30 GeV 3 / 2 and 
F + = 0.46 GeV 3 / 2 , corresponding to the results at zero order in the strong coupling a s . It should 
be noticed that for some of these parameters, as well as for the Isgur-Wise form factors £ and Tu2, 
the 0(a s ) corrections have been computed [|l5|, 17, 18|, 16]. However, since such corrections are not 



known for all the parameters needed in the present calculation, for consistency we use the values 
obtained at zero order in a s , including the effects of the known radiative corrections in the estimate 
of the theoretical uncertainty of our results. 

The universal form factors £, T\/ 2 and r 3 / 2 can be parametrized as follows: 

CM = (t^-) 2 (53) 
1 + U) 

t 1/2 (u) = 0.3 [1-0.5(^-1)] (54) 

r 3/2 (u/) = 0.3 [1 - 0.8(w - 1)] . (55) 



Eq.(53) is a useful parameterization of the Isgur-Wise form factor, obeying the normalization 
condition £(1) = 1 dictated by the heavy quark symmetry, and having a slope compatible with 
experimental data. The parameterizations for the Tj can be found, e.g., in [Q, |2(]], taking into 
account an uncertainty of 15% for the value at the zero recoil point oj = 1. pll 



Besides the already mentioned values for the strong coupling constants we use the following 
numerical values for the physical parameters appearing in the previous formulae: tub* —mB = 0.045 
GeV and m D *-m D = 0.142 GeV, V cb = 0.04 and V cd = 0.22 [|; a ~ 1.2 and c 2 = -0.2 §. As for 
the mass difference (5m 2 = mp* — mo, we use 5m 2 = 500 MeV [|J and the same value for <5m 2 = 
m Bl — m B- For the other quantities we use the theoretical determinations 5m = rriB — tub — 0.50 
GeV, T b = 0.30 GeV, T d = 0.14 GeV, T D » = 35 KeV. These values for the D* , D Q and B widths 
are consistent with the values for the strong coupling constants g and h given above (for a discussion 



see Ref. 

The decay width is given by 

r{rn B -m D ) 2 ps™ 11 * 

T(B° -> D + D-tt°) = ds„ ds+ , , (56) 

Jimo+rn^) 2 J s™ in ds + ds^ 
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where 

^' s '- BtD -* , ' = pi?i4 W <57) 

and 

9 

A=Y.Ak ■ (58) 

k=l 

Using the above parameters, and the formulae reported in the previous Sections, we find 

T{B° -> D + D-it°) ~ 5 x KT 16 GeV (59) 

and the corresponding branching ratio B(B° — ► D + D~n°) ~ 1 x 10~ 3 . Although it is difficult to 
assess the theoretical uncertainty related to this result, our calculation suggests that a relevant 
signal of the Cabibbo-suppressed B — > D + D~ir° decay should be detected at the B-factories. [] 

As for the size of the various terms appearing in the amplitude A, the 5*-type resonances in 
fig. la give a negligible contribution to the final result, whereas the contribution of the charmed 
intermediate states is dominant. Regarding the equal-time contribution, by itself this term would 
give a width of about 7 x 10" 17 GeV for process (ffl). Thus, although not quite dominant, it 
represents a contribution to the decay rate whose size is comparable to that of the main resonant 
terms. Moreover, its presence implies the specific dependence of the amplitude on s+ and s_ that 
takes the chiral symmetry constraint into account. 



We now consider the time-dependent decay probabilities {JSJ) and fllCD, which are given in terms 
of the functions Go, G c and G s in Eq.(||). Using the amplitudes A derived in Section ||| and the values 
of the input parameters listed above, we get the functions depicted in fig. 2. On the other hand, 
the functions obtained by the parameterization of the amplitude using the effective Lagrangian in 
Section || are depicted in fig. 3; the comparison between the two figures shows some agreement 
between the two methods of extrapolating the polar terms. 

We observe from both figs. 2 and 3 that the signal arising for the D* and Dq poles cannot 
be easily disentangled; concerning the contribution of the 2 + pole, the D% , we find that it is 
numerically small with respect to the other resonances, due to the small value of the form factor 
T 3 / 2 as compared to the constituent quark model value used in Ref. [Q]. 

For an assessment of the relative size of the various contributions to the time-dependent prob- 
abilities (||) and dl0|), we integrate the functions Go, G c , and ReG s and ImG s over the region 
bounded by Eqs.(^), ([?]), with s + > s_, corresponding to one half of the Dalitz plot. We find, using 
the functions in fig. 3: / G (s+, s^)ds + ds_ ~ 7 x 1(T 10 GeV 4 , / G c (s + , s^)ds + ds^ ~ -7 x KT 11 
GeV 4 , jReG s (s+, s_)ds + ds_ ~ 2 x KT 10 GeV 4 , / ImG s (s+, s^)ds + ds- ~ -8 x l(T n GeV 4 . 
These numbers may represent an indication on the sensitivity (hence, on the required statistics) 
of a combined time-dependent and Dalitz plot analysis of the events to the terms cos(Ami) and 
sin(Amt), as well as to sin(2/3) and cos (2/?). In particular, they suggest that the contribution of 
cos(2/3) to the time-dependent CP asymmetry may be sufficiently large to be identified. 



^Another Cabibbo suppressed B decay to charm mesons has been recently observed by the CLEO II Collaboration; 
it is the process B -» D* + D*~, with a measured branching ratio B(B -> D* + D*~) = [6.2±|;g ± 1.0] x 10~ 4 . H 
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Figure 2: The functions Go (up-left), G c (up-right), ReG s (down-left) and ImG s (down-right) in 
Eq.@ using the amplitudes reported in Section ||. The variables s± are in GeV 2 . 

6 The decay B -> D+D K s 

The kinematics of this process is quite similar to the case of the Cabibbo suppressed B — ► D + D~ir 
transition discussed above. Apart from the change m n — > mjf which reduces the available phase 
space, and the replacement of the D* resonance by the D* one, there are two important differences 
with respect to the case of the neutral pion in the final state. Indeed, as being induced at the quark 
level by the b — ► ccs transition, the amplitude for this channel is both color-allowed and Cabibbo 
favored by a factor Vcs/Vd, so that one expects an enhanced rate of events (and, possibly, a better 
efficiency of the Ks reconstruction as compared to the tt°). | Moreover, due to the flavour structure, 
beauty intermediate states (hence the amplitudes ^2,3,9 corresponding to fig. la) are absent, and 
the same is true for the D*~ resonances in Bq decays (amplitudes ^5.7 in fig. lb). In addition, the 
D* 2 does not contribute in the factorization approximation. Therefore, the resonant structure of 
*In this regard, also the uncertainty due to the penguin contributions should be reduced. 
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Figure 3: The functions Go (up- left), G c (up-right), ReG s (down- left) and ImG s (down-right) in 
Eq. @ using the parameterization of the amplitudes in Section El 



the amplitude turns out to be much simpler (although the parameters of the relevant Breit-Wigner 
forms are still to be measured). On the other hand, due to the significantly larger value of mx, the 
uncertainty implicit in the application of the chiral symmetry approach is expected to be larger 
for B — ► D + D~Ks- As far as SU(3)f breaking effects are concerned, we partially take them into 
account by replacing /„- — > fx (fx = 160 MeV) in the relevant formulae reported in the previous 
Sections. 

As a result, we obtain for B° -> D + D-R s the width r(2?° -> D + D-R s ) ~4x 1(T 15 GeV, 
corresponding to the branching fraction B(B° — > D + D~ K$) — 9 x 10~ 3 . This result indicates an 
enhancement of a factor 10 with respect to B — > D + D^ir°, rather than 20 naively represented by 
\V cs /V C( i\ 2 , and this is a consequence of the larger value of mx reducing the phase-space, and of the 
smaller number of intermediate resonances active in this case. 



The functions relevant for the time-dependent processes Eqs. (|[) and (10) are depicted in fig. 4 
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which shows the expected simpler Dalitz-plot structure, as far as the s + and s_ dependence is 
concerned, with respect to the process (|l]). 

The corresponding integrals over half of the Dalitz plot of such functions turn out to be: 
J G (s + ,s-)ds + ds- ~ 5 x 1(T 9 GeV 4 , / G c {s+, s^)ds + ds^ ~ -2 x 1(T 9 GeV 4 , 
J ReG s {s + ,s^)ds + ds- ~ 2 x lO" 9 GeV 4 , J ImG s (s + , s_)ds + ds_ ~ -6 x lO" 10 GeV 4 . 




Figure 4: The functions Go (up- left), G c (up-right), ReG s (down- left) and ImG s (down-right) in 
Eq.(H) for the transition B -> D + D-R s . 



7 Conclusions 

We have analyzed the three-body decays B — ► D Jr D~ir° and — > D + D~ Ks using a resonance 
model taking into account the constraints of the chiral symmetry on the relevant transition matrix 
elements. This method introduces a non-resonant, contact term as well as specific behaviour of the 
Breit-Wigner residue. 
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The main conclusion of our model can be seen in figs. 2 and 3, which show that the coefficient 
of cos(2/9) in the time-dependent rate, namely ImG s , is significantly different from zero over a 
sufficiently large portion of the Dalitz plot. Therefore, in principle, apart from its specific interest 
as a test of the chiral expansion in the heavy-quark theory, this channel might be useful to resolve 
the ambiguity in the determination of the CKM phase (3 {(3 — * 5-/3). Indeed, once sin(2/3) 
is measured from, e.g., B — > J/ipKs, the required information should be given by a suitable 
combination according to Eq.(|8|) of the two Dalitz plot distributions in the lower row of figs. 2 or 3. 

Qualitatively, our conclusions about the model dependence of the polar representation and 
the and tt° reconstruction efficiency agree with Ref. Q. From the numerical point of view, the 
differences with respect to j|] are mainly due to the inclusion of the equal-time commutator and the 
parameterization of the resonances. Indeed, the starting point of our calculation is the possibility 
of using the effective chiral lagrangian formalism for heavy mesons, offered by the smallness of the 
phase space available to the it . Another source of difference is the use of a smaller T3/2 form factor 
(as resulting from QCD sum rules) which depresses the contribution of D\. 

Finally, we emphasize the obtained large branching ratio B(B — ► D + D~ Kg) ~ 9 x 10~ 3 which, 
together with the simpler Dalitz plot structure and the better K$ reconstruction efficiency, can 
make this channel rather appealing for experimental analyses. 
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